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ABSTRACT 

We study the statistical distribution of extinct radio pulsars - isolated neutron stars, 
that just crossed the death line. An important element that distinguishes our study 
from other works is a consistent allowance for the evolution of the angle between the 
magnetic and spin axes. We find the distribution of extinct radio pulsars in spin period 
for two models: the model with hindered particle escape from the neutron star surface 
and the model with free particle escape. The total number of extinct radio pulsars is 
shown to be much smaller than for the model in which the evolution of an inclination 
angle is disregarded. The decrease in total number of extinct radio pulsars arises from 
the fact, that account for evolution of the inclination angle allows the transition of a 
neutron star to the propeller stage appear at much smaller periods (P ~ 5 — 10 s) 
than it is commonly assumed. 
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1 INTRODUCTION 



Previously, we showed (|Beskin Eliseeval[2QQ3h that in the 
case of free particle escape from the surface of an extinct 
radio pulsar for small distances d ~ 100 — 200 pc typical of 
the nearest neutron stars, their gamma-ray emission could 
be detected by the GLAST and AGILE missions. Fig. [1] 
shows that our estimate of the intensity of the emission 
from extinct radio pulsars for free particle escape from the 
neutron star surface a grees with the estimate obtained by 
([Harding, Muslimov Zhang 2002) for radio-emitting pul- 
sars to the numerical factor. The point is that the high- 
frequency emission is associated mainly with primary parti- 
cles and, therefore, it can also take place when no secondary 
plasma is produced. 

The mis match between our estim ate and the one done 
by ([Harding. Muslimov Zhandl2002l ) is actually the grav- 
itational factor. This difference arises due to the fact that 
we consider particle acceleration to be effective only at some 
distance from the star surface, at w hich the gravitational ef- 
fects can be neglected (|Beskin Elis eeva 2003). Let us draw 
provide more detailed consideration to this problem. 

Here we consider the model with free particle escape 
from the neutron star surface. Within this model the charge 
density of the particles flowing away from the polar cap dif- 
fers only slightly from the Goldreich- Julian charge density. 
General relativity effects become significant near the neu- 
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tron star surface ([Muslimov TsvganlFlQQd : lBeskinlll990ll . 
Since no particles are produced, the model of lMestell(|l999l ). 
in which no secondary plasma is generated, is appropriate in 
this case. The other model (Arons 1981) cannot be realized 
for extinct radio pulsars because it requi res a rev erse flow of 
secondary particles (for more details, see lBeskin|[T999l Here 
a short explanation should be done. The Mestel's model pre- 
sumes that none of the magnetic field lines near the pulsar 
surface are preferred because of the effects of general rela- 
tivity. Therefore the particles cannot be accelerated in this 
region ([Beskinlll99Qh . In other words, the longitudinal elec- 
tric field that arises from the mismatch between the plasma 
charge density and the Goldreich charge density decelerates 
rather than accelerates plasma particles. 

Nevertheless, the absence of regular acceleration does 
not imply that the plasma cannot fill the polar regions. Only 
effective particle acceleration is impossible. On the other 
hand, even at small distances from the pulsar surface r > ro, 
where 
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the general relativity effects become negligible. At such dis- 
tances the difference between the plasma charge density and 
the Goldreich charge density results in particle acceleration, 
at least on the half of the polar cap where the magnetic 
field lines deviate from the spin axis (i.e., those magnetic 
field lines for which the angle between the magnetic field 
line and the spin axis increases with distance from the neu- 
tron star). Below, we assume that particle acceleration be- 
gins only after the radius r = ro is reached. At smaller radii 
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Spin-Down Luminosity (erg/s) 

Figure 1. Predicted and observed dependencies of the gamma- 
ray intensity for radio pulsars on the rate of rotational energy 
losses under the assumption of free particle escape from the neu- 
tron star surface. Red line presents extinct radio pulsars. 

plasma on open magnetic field lines rotates rigidly with a 
neutron star. 

Now let us come back to the main point of this paper 
- the statistics of extinct radio pulsars. Since the intensity 
of the emission strongly depends on the spin period of the 
neutron star, radio-quiet pulsars with long spin periods (i.e., 
those that have long ceased to radiate in the radio frequency 
range or those with comparatively long spin periods at the 
stage of radio emission) can no longer be detected even under 
the assumption of free particle escape from the neutron star 
surface. Therefore, we can formulate three necessary condi- 
tions under which extinct radio pulsars are detectable: the 
validity of the model with free particle escape from the neu- 
tron star surface, a relatively small distance to such stars, 
and short spin period of a neutron star. 

Clearly, the observability of extinct radio pulsars de- 
pends on their total number in the Galaxy. Knowing the 
total number of extinct radio pulsars and assuming the 
equiprobable distribution of these stars in the Galaxy, we 
can estimate the number of extinct radio pulsars can be 
within about 200 pc of the Sun. Another no less important 
parameter is the spin period of a neutron star at which a 
star passes to the stage of an extinct radio pulsar. 

In this paper we consider the stationary distribution 
of extinct radio pulsars within the scope of two models: 
the model w ith hindered particle escape fr om the neutron 
star surface ([Ruderman fc Sutherlandll 19751 ) and the model 
with free particle escape (|Aronslll979l ). An important ele- 
ment that distinguishes our study from other results is a 
consistent allowance for the evolution of an angle between 
the magnetic axis and the spin axis. Throughout the paper 
we assume that no neutron stars are born at this stage. Fi- 
nally, we disregard the possibility of the magnetic field decay 
during the stages of active and extinct radio pulsars. 

As a result, we find the distribution of extinct radio pul- 
sars in spin period. We show that total number of extinct 
radio pulsars is much smaller than the number of extinct 
radio pulsars within the scope of model without the evolu- 
tion of the inclination angle. The reason for such a decrease 
in the total number of these objects is the following: if the 



evolution of an inclination angle is taken into account, the 
transition to the propeller stage occurs at much shorter spin 
periods (P ^ 5 — 10 s) than it is possible for the standard 
model (P ^ 100 s), for which the evolution of an inclination 
angle is commonly disregarded. 



2 BASIC EQUATIONS 

Our goal is to determine the number of neutron stars that 
are no longer radio pulsars, but have not yet passed to the 
propeller stage. In other words, we seek the distribution of 
neutron stars in the interval Pcz(x,P) < P < Pprix^B)^ 
where the period Pdix^B) corresponds to the death line 
of radio pulsars, and the period Ppr{x,B) corresponds to 
the transition of neutron stars to the propeller stage (see 
Fig. [2] and [3)). We take into account the fact that both 
P = Pd(x,P) and P = Ppr{x,B) depend on the angle of 
axial inclination x (this is a new element that has not been 
studied previously). 

The death line of radio pulsars can be written as follows 
(jBeskin et al.lll993h 

P4x)-BT{cosxr-''s. (2) 

It is known that the transition to the propeller stage oc- 
curs when the Schwarzman radius (i.e., the radius at which 
the pressure of the ambient medium is equal to the mag- 
netodipole radiatio n pressure) is t he sam e as the radius of 
the Alfven surface iLipunov et al.1 (|l996h . Therefore, to de- 
termine the x-dependence of the transition period Ppr we 
note that the magnetic dipole losses 

T T OO 1 BqQ^R^ . 2 /Q\ 

Lmd = -Jr^il = ^ sm X (3) 
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(Po is the magnetic field at the pole of a neutron star with 
a radius R and a moment of inertia Jr) are proportional to 
the combination Bosinx- As a result, we obtain 

Pprix) = Pe sin^/^ X, (4) 



wher e the expressi on for Pe has the standard form 
(Lipun ov et al.lll996h 




Here, poo and Coo are the density and the sound speed in the 
interstellar medium, respectively (cr is in units of 10^ cm 
s~^, Pext — ^TTpoocio IS the corrcspoudiug energy density 
of the magnetic field {{Bext)-6 is in units of 10~^ G). We 
also assume that the velocity of neutron stars Voo is much 

G reater th an the sound speed in the interstellar medium 
Popov et al. 2000). 

To determine the stationary distribution function of ex- 
tinct radio pulsars, we introduce their distribution function 
N{P, X, B) in spin period P, inclination angle x? and mag- 
netic field B. The kinetic e quation can be presented in gen- 
eral form as (Bes kin et al.lll99 3) 

dP\St )^ dx\5t ) ^ ' 
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Figure 2. Isolated neutron star evolution under the assumption 
of hindered particle escape. Regions I and II correspond to active 
and extinct radio pulsars, respectively. The death line of active 
radio pulsars is shown for three different magnetic field strengths. 
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Figure 3. Isolated neutron star evolution under the assumption 
of free particle escape. Regions I and II correspond to active and 
extinct radio pulsars, respectively. 



In other words, we ignore the magnetic field decay over 
the lifetime of an extinct radio pulsar the same way as we 
did for radio pulsars. The validity of the assumption for 
normal radio pulsars seems proven, since their lifetime does 
not exceed the characteris tic magnetic-field evolution time 
(jPopov h Prokhorovllioooh . As we show below, if the evo- 
lution of the angle of inclination is taken into account, the 
lifetime of a neutron star during the extinct radio pulsar 
stage is also short enough. 

The derivatives 5P/5t and Sx/^t (|6|) are known func- 
tions of the magnetic field B, the period P, and the an- 
gle of inclination x- Further, if hindered particle escape 
takes place and the spin period of a pulsar P > Pd{x^B), 
plasma fills only the inner region of the pulsar magneto- 
sphere (Isto min M osvagin 1995). Thus, for the model with 
hindered particle escape from the radio-pulsar surface, we 
assume that the energy losses closely match the magnetic 
dipole losses (|3|): 
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St ^ 
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-10 



'^^B^oP- 



1 . 2 

sm X, 



^ sin X cos X s ^ . 



(7) 

(8) 



It should be emphasized here that within the scope of 
hindered particle escape model the following ratio is invari- 
ant for each individual neutron star throughout the stage of 
extinct radio pulsars. 



id = 



cos X 



(9) 



As a result, an inclination angle of each extinct radio pul- 
sar decreases during the evolution (providing the validity of 
hindered particle escape model). Moreover, an inclination 
angle chi decreases at the same rate as the spin period of 
the neutron star increases. 

As for free particle escape model, here we assume that 



the neutron star magnetosphere is still completely filled with 
plasma that screens the longitudinal electric field. In this 
case, the magnetic dipole losses are known to be completely 
screened. Therefore, all the energy losses must be associ- 
ated with the long i tudinal curre nts in the magnetosphere 
(|Beskin et al.lll993( . iMestell Il999l) . In other words, the en- 
ergy losses are the same as those for radio pulsars. Thus, 
within the sc ope of free particle escape model we obtain 
(iBeskin et al.| [l993) 



1 n-15 d10/7 7-,l/14 2d 
10 P COS X, 



^ ^ 10-^^B,T^^/''cos^^-^X 



(10) 
(11) 



Here 2(i 1.5 is a dimensionless parameter. So, for the 
model with free particle escape the ratio 



smx 



(12) 



is invariant throughout the time a neutron star spends as an 
extinct radio pulsar. As we can see, both an inclination angle 
X and a spin period increase at the same rate providing the 
particles freely escape a neutron star surface. An inclination 
angle for this model tends to 90°. So, for both models of 
particle escape function the angles tend to pass into the 
region where the energy losses are minimumal. 

Further, function U on the right-hand side of Eq. ^ 
is the source of extinct radio pulsars. To make a reason- 
able assumption about this function we should recall how it 
is defined in case of active radi o pulsars . Usin ^ the results 
obtained in the monog raph iBeskin et al.1 (|l993l ) , we assume 
that the birth function of neutron stars can be presented as 
a product of the probabilities of their initial distributions in 
angle Xi period P, and magnetic field B: 



U ^U{x)Up{P)Ub{B). 



(13) 



Unfortunately, we have no information about the dis- 
tribution of nascent neutron stars in axial inclination x C^)- 
Therefore, we have to assume that the distribution of neu- 
tron stars in their initial inclination angles x is equiprobable, 
so that 

2 



U{x)^ 



(14) 



In order to determine the birth function U{P), we find 
the period distribution for pulsars with characteristic ages 
td = P/2P greater than the kinetic ages Tkin = z/v±, which 
can be determined for pulsars with known proper motions. 
Here v± is the pulsar velocity perpendicular to the Galac- 
tic plane. We use the fact that the thickness of the disc 
of radio pulsars, Zp ^ 300 pc, is much larger than the 
thickness of the disc o f supe rnova remnants, Zsn ~ 10 pc 
iMarochnik Suchkovl (|l996l ). where most of the neutron 
stars are assumed to be born. As a result, the radio pul- 
sars from this sample have not changed their periods signif- 
icantly. Therefore, the distribution of such pulsars in spin 
period P should be close to the true birth function of neu- 
tron stars U{P). 

Fig. |4] generally matches the distribution obtained pre- 
viously (jBeskin et al.lll993l ). The distribution seems to be 
nearly equiprobable. At least, it has no distinct maximum 
at short periods. Therefore, in case of active radio pulsars 
we can assume that 
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Figure 5. Distribution of young radio pulsars {Q < 1) in mag- 
netic field Bi2 on a logarithmic scale. Data from the ATNF cat- 
alogue are used to construct the distribution. 



Figure 4. Spin period distribution for 37 active radio pulsars 
with r£, > Tkin- Data from the ATNF catalogue are used to 
obtain the distribution. 



Up{P) =:Up = const. 



(15) 



We use this distribution function to analyze the statis- 
tical distribution of radio pulsars. On the other hand, we 
see that a sharp (almost to zero) decrease in the number of 
radio pulsars is observed at periods P > 1 s. At the same 
time, the decrease in the number of radio pulsars for the 
complete sample in the period interval 1 < P < 2 s is not 
so essential. Therefore, we assume that no neutron stars are 
born at the stage of extinct radio pulsars: [7 = 0. 

Finally, magnetic-field distribution function of radio 
pulsars can be written as (?) 

t/.(B) = ^;j^(^(|)^xp(-S/Bo), (16) 

where Po = 10^^ G and 7 = 2. In order to refine the mag- 
netic field dependence of the distribution function of radio 
pulsars, we use the logarithmic chart plotting the number of 
radio pulsars N{B > B p) vs . magnetic field P for 1348 pul- 
sars (jManchester et al.ll2005l ). The pulsars are divided into 
two groups by the dimensionless parameter Q that reflects 
their basic characteristics: 

/ P \ 11/10 / P 



The pulsars with Q > 1 and Q < 1 lie near a nd far from 
the death line (for more details, see lBeskin et al ](l884). Fig- 
ures [5] and [6] show the distribution functions of extinct radio 
pulsars with Q < 1 and Q > 1, respectively. Surprisingly, 
despite the fact that the number of observed radio pulsars 
increased more than 3 times over the last ten years, we ob- 
tained almost the same results that were found previously 
for the statistics of about 400 pulsars. As a result, the distri- 
bution fu nction of radio pulsa rs with Q < 1 can be written 
as follows iBeskin et all (|l993h 



: 2kNNfP 
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(18) 



sin^ X 



'x X e. 



Here, fciv « 4.4 is a normalizing factor, Q = 9[Pi{x) — P], 
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Figure 6. Distribution of old radio pulsars {Q > 1) in magnetic 
field Bi2 on a logarithmic scale. Data from the ATNF catalog are 
used to construct the distribution. 



and Nf is a number of active radio pulsars far from the 
death line. 

Thus, we consider the stationary distribution of extinct 
radio pulsars without neutron stars being born at this stage 
and with the magnetic- field evolution disregarded. Accord- 
ingly, the kinetic equation (|6]) can be rewritten as 

The solution of the kinetic equation for the model with hin- 
dered particle escape is given by 



N{P,x) 



(20) 



sinx 

For the model with free particle escape from the pulsar sur- 
face it is 

F(L) 



NiP,x) 



COSX 



(21) 



To solve the kinetic equation for the two models, we should 
define the boundaries of the region of extinct radio pulsars. 
For the model with hindered particle escape from the neu- 
tron star surface (the Ruderman- Sutherland model) these 
boundaries are determined by the death line of radio pulsars 
(Gi = ^[Pi(P, x) ~ P]) and by the line where the transition 
of extinct radio pulsars to the propeller stage takes place 
(G2 = 0[sinx — P'^/Pe])- In within the scope of free parti- 
cle escape model, the transition to the region in which an 
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inclination angle is close to 90 ° (jBeskin Nokhrinal[2003 l 
should be added to the above boundaries. 

Let us designate the regions of active and extinct radio 
pulsars as regions I and II, respectively. To determine the 
function F{I), we should match the solutions at the bound- 
ary between the regions of active (region I) and extinct (re- 
gion II) radio pulsars. Due to the fact that the neutron star 
fluxes on both sides of the death line should remain equal, 
we obtain for the model with hindered particle escape 

. / dxi dh/dP \ 
' I dP dic/dx ) 

' V dP did/dx ) 

Here the period and angle derivatives correspond to the 
death Une of pulsars P = Pd{x)- At the same time, for the 
model with free particle escape, we obtain 

N2[Pd{x),x] = Ni[Pd{x),x]- (23) 

As a result, in the Ruderman-Sutherland model, we can 
approximate the function F(I) by the following expression 




Figure 7. Distribution of extinct radio pulsars N2(P) in the 
scope of hindered particle escape model with (solid lines) and 
without (dashed line) allowance for the evolution of the inclina- 
tion angle. 



F{Ia) « AiBi2 



cos X 



(24) 



where n ^ -0.6, and A{Bi2) oc ^12(1 + B12) 



3 RUDERMAN-SUTHERLAND MODEL 

Let us first consider the statistics of extinct radio pulsars 
assuming the model of hindered particle escape from the 
neutron star surface is valid. As it was mentioned above, 
the particle work function is fairly large in this case, and 
particles fill only the equatorial regions of the magnetosphere 
with the strongly curved magnetic field, in which secondary 
particles can still be produced. Therefore, the energy losses 
for the Ruderman-Sutherland model can be described with 
a sufficient accuracy by a magnetic dipole mechanism. 

Let us now turn to Eq. (|22|) derived from the conserva- 
tion condition for the neutron star flux on the death line of 
radio pulsars. This expression deflnes the function F(I) and, 
hence, the distribution of extinct radio pulsars in period P, 
magnetic fleld B, and inclination angle x for the model in 
question. In this equation, the derivative Pi corresponds to 
the region of active radio pulsars and can be specifled by 
iBeskin et af ] ( jlQQSl ): 



SP^ 

St 



10' 



-15^10/7pl/14 



1.5 

COS X- 



(25) 



At the same time, the derivative P2 is specifled in the 
region of extinct radio pulsars, and Eq. JT]) should be used 
here. Next, let us recall that the region of extinct radio pul- 
sars in the scope of hindered particle escape model is limited 
by the death fine of active radio pulsars (|2]) and the line of 
the transition to the propeller stage (|4]). 

As a result of the above expressions and using the fact 
that distribution function for normal radio pulsars is deflned 
by Eq. {[19]), Eq. (|22)) can be rewritten as 



N2[Pd{x),x] ^ 0.7kNNfB%\l + B12 

1 _ ^^o^/ 



-3.7 



(26) 



COS 



sin^X 



x(3.5 + t/x)- 



On the other hand, the solution of the kinetic equation 
for the model with hindered particle escape providing the 
distribution of extinct radio pulsars is stationary, no stars 
are born at this stage, and the magnet ic-fleld evolution is 
disregarded, can be represented by Eq. (|20p . Thus, for the 
distribution of extinct radio pulsars on the death line, we 
obtain 



N2[Pd{x),x] 



F(cos^-'\/bT) 



(27) 



smx smx 

Further, we rewrite the distribution function for the ex- 
tinct radio pulsars derived from the conservation condition 
for the flux through the death in the same form as Eq. 1271 
For this purpose, we use the following change of variable: 



B 



,8/15 



(28) 



As a result, the distribution function for the extinct 
radio pulsars can be flnally transformed to 

iV2[P,B,x] = 0.7A:ivA^/Bi\/^(l + Bi2)-'-^ x (29) 
^G{cosx/P), 



smx 



-6162, 



where 



G{0 



:i<i)(3.5 + 4®)c^"(6. 



Here we use two auxiliary functions: 



C(0 = ^12 ^ 



siO = [^-ciOr' = sinx 



cosx, 

1/2 , 



(30) 

(31) 
(32) 



Figure [71 shows the distribution of extinct radio pulsars 
in spin period P under the assumption of hindered particle 
escape for three different values of limiting spin period Pe 
([6|). The results of numerical calculations shown in this flg- 
ure are in good agreement with the analytical estimate (|24p 
obtained for the Ruderman-Sutherland model. As we see, if 
the evolution of an angle of inclination of the magnetic axis 
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to the spin axis is consistently taken into account, total num- 
ber of neutron stars at the extinct radio pulsar stage proves 
to be much smaller than the same total number within the 
model for which the evolution of an inclination angle is dis- 
regarded. On the other hand, the number of pulsars with 
comparatively short periods P ^ 2 — 4 s at this stage is 
even larger than that for the standard model (the dashed 
line in Fig. [7|). This is attributable to the fast transition of 
neutron stars to the region of small angles of inclination, 
where they are not only accumulated through the decrease 
in spin-down rate, but also disappear due to the transition 
to the propeller stage. 

Finally, it is very important that, if the evolution of the 
angle of axial inclination is taken into account, the transi- 
tion to the propeller stage can also occur at short neutron 
star spin periods P ~ 5 — 10 s, or even at P ~ 1 — 2 s (see 
Fig. [2}. Indeed, as we mentioned above, the transition to the 
propeller stage occurs when the Schwarzman radius, which 
can be determined from the balance between the magne- 
todipole and external radiation pressures, becomes equal to 
the Bondi-Hoyle radius. However, since the magnetodipole 
losses at small angles of axial inclination are essentially sup- 
pressed, the transition to the propeller stage takes place even 
at short neutron star spin periods: 



Psisinx) 



1/2 



(33) 



which is about 5 — 10 s at Pe ^ 100 s (|6|). Recall once again 
that if the model with hindered particle escape from the 
neutron star surface is valid, extinct radio pulsars cannot be 
detected with currently available instruments. 



4 ARONS MODEL 

Let us now discuss the statistics of extinct radio pulsars 
assuming the model with free particle escape from the neu- 
tron star surface is valid. Recall that for this model, the 
energy losses for both extinct and active radio pulsars are 
the same and are associated with the longitudinal currents 
in the magnet osphere. Wa want to stress here that within 
the model of free particle escape gamma-ray emission from 
extinct radio pulsars with short (several seconds) spin pe- 
riods and providing relatively small distances to such stars 
can be detected. 

Let us consider Eq. (|23)) . which is derived from the con- 
dition for the equality between the neutron star fluxes on 
both sides of the death line. Using the distribution function 
of active radio pulsars ([T9|) derived previously Beskin ~e't al.l 
((1993) and refined by using the new data in this paper, we 
obtain the following expression for the distribution function 
of extinct radio pulsars on the death line: 



A^2[Pd(x),x] = 0.6A:iviV/PiV(l + Pi2) 



(34) 



COSX -0.21 

— COS X- 



On the other hand, the solution of the kinetic equation 
for the model of free particle escape from the neutron star 
surface can be represented as (|21|) . Accordingly, for the pul- 
sars on the death line, the solution of the kinetic equation 
is given by 



N2[P4x),x] 



F{U) _ F[sinx/(B%''cos°-'^x)] 



(35) 



COSX COSX 
Recall that, as for the Ruderman-Sutherland model, 
this solution is obtained assuming the distribution of ex- 
tinct radio pulsars to be stationary. At the same time, we 
disregarded the possibility of stars being born at this stage 
and the possibility of magnetic field decay. 

Next, it is necessary to consider the boundaries of the 
region of extinct radio pulsars for the model with free parti- 
cle escape. This region is limited by the death line of radio 
pulsars, the line of the transition of pulsars to the propeller 
stage, and the line of the transition of pulsars to the region 
where the angle of inclination of the magnetic axis to the 
spin axis of a neutron star is close to 90° (see Fig. [3]). As in 
the case of hindered particle escape, the death line of radio 
pulsars is given by Eq.©. The transition of extinct radio 
pulsars to the propeller stage is defined as follows: 

Pprix) = Pecos^^\. (36) 

Here, the expression for the limiting period Pe also 
has the standard form (|6|). Finally, the third boundary of 
the region of extinct radio pulsars is the line of the transi- 
tion to the region where an inclination angle approaches 
90°. For this boundary we have the following expression 
(jBeskin Nokhrina.2004 ): 



COSX 



2^i?\V2 



(37) 



Unfortunately, within the scope of free particle escape 
model, the distribution function of extinct radio pulsars de- 
rived from the flux conservation on the death line cannot be 
analytically rewritten in the same form as ()35|) . Still, we can 
find analytical approximations for the large or small incli- 
nation angles. Thus, for small angles of axial inclination, 

1 



P(/) = -A:iviV/Pi2 (1 + ^12) 



(38) 



At the same time, for angles of inclination close to 90°, 
we obtain the following expression: 



F{I) = -kNNfB^^2\l + Bi2) 



-3.7j-1.7 



(39) 



Now, using the expressions obtained for large and small 
angles of inclination x, we can approximate the function 
F{I) for all inclination angles < 0% < 7r/2 as 

0.51 



F{I) 



(40) 



(l + 0.1Pl-7pl)0-9' 

Finally, the distribution of extinct radio pulsars for the 
model with free particle escape can be represented by the 
function 



N2[P,B,x] = 0.3kNNf X 

i^lV(l + i^l2)-^-' 



(41) 



rOi 



cosxil + 0.1Pi-7(smx/P)3-i)0-9 ' 

Recall that the normalizing factor in this expression is 
/cat ~ 4.4, and Nf gives the number of normal radio pulsars 
far from the death line. 

Figure [8] shows the distribution of extinct radio pulsars 
in spin period P under the assumption of free particle escape 
from the neutron star surface. As we can see, for this model, 
just as for the Ruderman-Sutherland model, consistent al- 
lowance for the evolution of the angle between the magnetic 
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Figure 8. Distribution of extinct radio pulsars N2{P) within the 
scope of free particle escape model. 

axis and the spin axis causes the total number of extinct 
radio pulsars to decrease. Such a decrease appears because 
neutron stars can pass to the propeller stage at much shorter 
spin periods than assumed previously. Indeed, if inclination 
angles are close to 90°, extinct radio pulsars can pass to the 
propeller stage even at periods 

P^^^P^^cosx)^^^. (42) 

Thus, if an angle of inclination of the magnetic axis to 
the spin axis is close to 90°, the spin period of a pulsar 
passing to the propeller stage could be P ^ 5 — 10 s. 



5 CONCLUSIONS 

We found the theoretical distribution of extinct radio pul- 
sars in spin period within the scope of two models: the model 
with hindered particle escape from the neutron star surface 
and the model with free particle escape. In both cases, we 
consistently took into account the evolution of an angle be- 
tween the magnetic axis and the spin axis. 

As a result, we showed that for both models all the 
details of the calculations, such as the death line of active 
radio pulsars, the transition of extinct radio pulsars to the 
propeller stage, and the evolution of extinct radio pulsars, 
actually depend severely on the evolution of an inclin at ion 
angle . Unfortunately, with a few exception s ([Beskin et al.l 
1 1991 iRedmbau &: de Freitas Pachecol I2OO1I ) . no consistent 
analysis of the evolution of an inclination angle has been 
performed previously. 

As we proved, such an allowance for the evolution of 
the axial inclination in analyzing the statistics of extinct 
radio pulsars causes the total number of such radio pulsars 
to decrease, thereby making it difficult to detect them with 
currently available instrumentation (the GLAST or AGILE 
missions). The decrease is due to the possibility for extinct 
radio pulsars to pass to the propeller stage at shorter spin 
periods P ~ 5 — 10 s than it is possible for the model without 
account of inclination angle's evolution. Nevertheless, if we 
consider the model with hindered particle escape and take 



into account the evolution of an angle the number of 
extinct radio pulsars with particular short spin periods (P ^ 
2 — 4) s is even larger than that in the standard model. 
Finally, it should also be stressed here that the detection 
of extinct radio pulsars would be direct evidence for free 
particle escape from the neutron star surface. 



6 ACKNOWLEDGEMENTS 

We are grateful to S.B.Popov and M.E.Prokhorov for fruitful 
discussions and numerous advice. This work was supported 
by the Russian Foundation for Basic Research (project nos. 
02-02-16762 and 05-02-17700). 



REFERENCES 

Arons J., 1979, Space Sci. Rev., 24, 437 

Arons J., 1981, Astrophys. J., 248, 1099 

Beskin V.S, 1990, Pisma Astron. Zh., 16, 665 [Sov. Astron. 
Lett. 16, 286 (1990)] 

Beskin V.S, 1999, Phys. Uspekhi 40, 659 

V.S. Beskin and S.A. Eliseeva, 2003, Astron. Lett., 29, 20 

Beskin V.S., Gurevich A.V., Istomin Ya.N., 1984, Astro- 
phys. Space Sci., 102, 301 

Beskin V.S., Gurevich A.V., Istomin Ya.N., 1993, "Physics 
of the pulsar magnetosphere" , Cambridge University 
Press 

V.S. Beskin and E.E. Nokhrina, 2004, Astron. Lett., 30, 
685 

A. K. Harding, A. G. Muslimov and B. Zhang, 2002, As- 
trophys. J., 576, 366 

Ya. N. Istomin and D. V. Mosyagin, 1995, Astron. Zh., 72, 
826 [Astron. Rep. 39, 735 (1995)] 

V.M. Lipunov, 1987 "Astrophysics of Neutron Stars", 
Nauka, Moscow 

V.M. Lipunov, K.A. Postnov and M.E. Prokhorov, 1996, 
Astrophys. Space Phys., 9, 1 

Manchester R.N., Hobbs G.B., Teoh A., Hobbs M., 2005, 
AJ, 129, 1993 

Marochnik L.S., Suchkov A.A., 1996, "The Milky Way 

Galaxy", Gordon and Breach 
Mestel L., 1999, "Stellar Magnetism", Oxford University 

Press 

A.G. Muslimov and A.I. Tsygan, Astron. Zh., 1990, 67, 263 

[Sov. Astron. 34, 133 (1990)]. 
S.B. Popov, M. Colpi and A. Treves, 2000, Astrophys. J., 

530, 896 

S.B. Popov and M.E. Prokhorov, 2000, Astron. Astrophys., 
357, 164 

Regimbau T., de Freitas Pacheco J. A., 2001, A&A, 374, 
182 

Ruderman M.A., Sutherland P.G., 1975, ApJ, 196, 51 
Tauris T.M., Manchester R.N., 1998, MNRAS, 298, 625 



